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Abstract 

In models with flavor symmetries in the leptonic sector leptogenesis can take 
place in a very different way compared to the standard leptogenesis scenario. We 
study the generation of a B — L asymmetry in these kind of models in the flavor sym- 
metric phase pointing out that successful leptogenesis requires (i) the right-handed 
neutrinos to lie in different irreducible representations of the flavor group; (ii) the 
flavons to be lighter at least that one of the right-handed neutrino representations. 
When these conditions are satisfied leptogenesis proceeds due to new contributions 
to the CP violating asymmetry and — depending on the specific model — in several 
stages. We demonstrate the validity of these arguments by studying in detail the 
generation of the B — L asymmetry in a scenario of a concrete A^ flavor model 
realization. 
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1 Motivation 



Observational data from the abundances of light elements (D, 3 He, 4 He and Li) in addition 
to precision observations of the cosmic microwave background (CMB) temperature fluc- 
tuations allow the determination of the cosmic baryon asymmetry, Y& B = (rig — n^)/s = 
(8.75 ± 0.23) x 10~ n (with tib (n§) the baryon (antibaryon) number density and s the 
entropy density) pp. Though the conditions for dynamically generating this quantity are 
well known and established |2j the cosmic baryon asymmetry poses a puzzle in particle 
physics: the standard model (SM) fails to explain such a large asymmetry, thus implying 
the presence of new physics accounting for Y^ B . 

Leptogenesis is a scenario in which a lepton asymmetry Y& L is generated in the lepton 
sector and partially reprocessed into Y& B by SM electroweak sphaleron processes (for a 
comprehensive review see [3]). The generation of Y"a l requires, in addition to CP violation 
and departure from thermodynamical equilibrium, lepton number breaking. Accordingly, 
in these class of scenarios two in principle unrelated puzzles are linked, the origin of 
neutrino masses and the baryon asymmetry. Among all the possible neutrino mass models 
present in the literature the standard seesaw (type I seesaw) |4J provides the framework for 
standard leptogenesis, in which the lepton asymmetry proceeds via the out-of-equilibrium 
and CP violating decays of the lightest right-handed (RH) electroweak singlet neutrino. 

Most of the studies of leptogenesis are based on the assumption that there is no new 
physics between the lepton number breaking scale and the electroweak scale that can 
sizable affect the way in which leptogenesis takes place. Though some analysis in scenarios 
including flavor symmetries above the electroweak scale have been done, and have proved 
that the presence of new energy scales and new degress of freedom may have an impact on 
the way leptogenesis proceeds [5j El [7J [8] , all of them are based on the same assumption, 
namely the lepton number breaking scale is below the scale at which the flavor symmetry 
is broken B 

The idea to ascribe to a flavor symmetry to explain particle masses and mixings dates 
backs to the late 1970's [UJ. Originally flavor symmetries were introduced to explain 
quark structures and only after neutrino oscillation data the use of horizontal symmetries 
1 The exceptions being references [9| flO] . 
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in the lepton sector has become more challenging and interesting. In particular, in the last 
years it has been shown that lepton mixing may be well described by discrete non Abelian 
symmetries (see [12J and references therein for further details). Given that in these kind 
of models both, the lepton number and flavor breaking scales are free parameters the 
question about how does leptogenesis proceeds in the flavor symmetric phase proves to be 
quite reasonable. 

In more detail, let us suppose to have a non Abelian flavor symmetry group Gp un- 
der which the RH neutrinos and left-handed SM leptons have definitive transformations. 
We introduce a cutoff scale A since we will deal with non-renomalizable operators. The 
scenario which we are interested in is the following: the lepton number breaking scale, 
characterised by the RH neutrino mass, Mjy, is larger than the scale at which Gp is bro- 
ken, vp. This means that the Yukawa mass matrices above and below vp are different: in 
particular the Yukawa Dirac matrix below vp is proportional to the Dirac mass matrix — 
the proportionality factor represented by the SM Higgs vacuum expectation value (vev). 
Flavon masses, M^, are taken as free parameters, clearly not too far from vp but above 
it. Clearly it holds A > M N , M^> v F . 

Our discussion is based on the class of symmetries that explain neutrino masses and 
lepton mixings and may be be generalized to any flavor symmetry, Abelian or not, discrete 
or continuos. We will exemplify our arguments by doing a full analysis of the generation 
of the B — L in a concrete model that at low scale exhibits exact TriBiMaximal (TBM) 
mixing at leading order. The reason for this choice is very simple: it has been shown 
that when exact TBM mixing is induced by type I seesaw the CP violating asymmetry 
is zero and acquires a non-vanishing value only when lepton mixing deviates from TBM 
[T3l [HI [7] . In type II seesaw models this could be also the case if only two electroweak 
triplets are present and they both are family singlets [16]. Extended models featuring 
interplay between type I and II seesaws have been also analysed and the conclusion is that 
in these cases it is also possible to identify a generic class of minimal models in which the 
CP asymmetry vanishes as well [13] . 

The most recent analysis indicates that TBM is not anymore in perfect agreement with 
neutrino experimental data since the reactor angle predicted by TBM is zero, while this 
value is now excluded at 3a level [IT]. However TBM remains a good approximation for 
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the lepton mixing matrix and we will consider a model that predicts exact TBM for its 
simplicity in showing the feasibility of leptogenesis in the regime A > Afjy, > Vp. 

The paper is organized in the following way: next section is general and we enumerate 
the general conditions necessary to obtain a CP asymmetry, ejv, in the flavor symmetric 
regimen. Section [3] shows how the proposed conditions work by mean of a complete 
analysis of a specific model based on the flavor symmetry The model main features and 
neutrino phenomenology are briefly discussed and the generation of the B — L asymmetry 
is explained in detail. Section [4] is devoted to our conclusions. The calculation of the 
reaction densities necessary for the analysis of the washout processes studied in sec. [3] are 
given in appendix |A"{ 

2 Leptogenesis in the flavor symmetric phase 

We have already anticipated in the introduction that in models for leptonic flavor mixing 
four energy scales can be distinguished, namely a cutoff scale A — or in general a scale 
of heavy matter — , the lepton number breaking scale — determined by the RH neutrino 
masses — M^, the flavons scale — determined by the masses of the scalar s that trigger 
the flavor symmetry breaking — and the scale at which the flavor symmetry is broken, 
Vp. Though A > Mn, M<p,VF the scales Mjy, and Vp, being free parameters, can 
follow any hierarchy. Since we are concerned about leptogenesis in the flavor symmetric 
phase it is clear Vp < Mjv, M^. This constraint in turn has an implication: if the flavor 
symmetry enforces the RH neutrinos to belong to the same Gp irreducible representation, 
R, leptogenesis will not be achievable: in the flavor symmetric phase the RH neutrinos 
have a common universal mass and therefore the CP violating asymmetry €n vanishes 

As a consequence, in these kind of models viable leptogenesis requires RH neutrinos 
to belong to different Gp irreducible representations Ri, so a mass splitting among the 
masses of the different representations can be accomodated. Let us assume the existence 
of k electroweak lepton doublets placed in r representations L r , m RH neutrinos lying 
in p representations N p and n electroweak singlet scalars arranged in q representations 
S q . Assuming the SM Higgs SU(2) doublet, H, to be a Gf singlet the z-th RH neutrino 
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representation can only decay to final states containing Lj. Accordingly, three type of 
models can be distinguished: 

1. For any RH neutrino representation iVj ~ Ri there is a lepton doublet representation 
Li ~ R[ with which a gauge flavor invariant renormalizable operator Li Ni H can be 
built 

2. Only for a set of the RH neutrino representations a gauge flavor invariant renormal- 
izable operator Li Ni H exist 

3. For non of the RH neutrino representations the operator Lj Ni H can be built. 

In cases [T] and [2] the standard one-loop vertex and wave-function corrections to the tree- 
level decay exist however the CP violating asymmetry derived from their interferences and 
the corresponding tree-level process vanishes. The proof of this statement is easy. Suppose 
Ni transforms as the Ri representation of Gp. To recover the correct kinetic term we know 
R*Ri = 8 ai p. should hold (here a», $ are the indices of the Ri representation). The same 
applies for Lj ~ R^. Now, the standard contribution to the CP asymmetry involves the 
imaginary part of the off-diagonal elements of the matrix YY\ where Y is a generic 
Yukawa coupling matrix defined in the stage of unbroken flavor symmetry. This Yukawa 
coupling matrix is given by the Clebsch-Gordan coefficients arising from the contractions 
R'iRi, thus implying that the matrix YY^ is determined in turn by the contractions 
R^RiR^R'*, which demonstrates that YY^ is diagonal. Viable leptogenesis is, therefore, 
possible only if new contributions (e^ ew ) to e^v are present, and this is possible only if at 
least for one of the RH neutrino representations the condition > is satisfied ^ 

2 Here i labels the index representation, not the flavor index and with 2Vj ~ Ri(Li ~ R^) we indicate 
that Ni(Li) transforms as the representation Ri(R' i ). We introduce the generic R[ representation for Li 
to be as general as possible, defining R\ as the representation that contracted with Ri has in its Clebsch- 
Gordan series a singlet of Gf- Indeed for the discrete groups with real triplet or doublet representation 
it holds that R^ — Ri. For these kind of groups it makes sense defining the R* only when the physical 

field-a scalar or a fermion-is complex. 

3 This condition, in addition with the mass splitting among the different RH neutrino representations, 

guarantee that the novel loop corrections to the tree-level decay contain an imaginary part, an essential 

requirement for e^ cw ^ 0. 
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In case [3] it is possible that iVj may have n > 2 body decays by means of non renormal- 
izable operators. In this scenario one should modify the standard case in order to include 
n-body decays. However, due to space suppression factors the CP asymmetry generated 
is expected to be small. For this reason in what follows we do not consider this case. 

Assuming the flavons are lepton conserving states no B— L asymmetry can be generated 
via decays. However, once the condition Mn > is satisfied they can play an essential 
role in the generation of the B — L asymmetry, not only because they lead to novel 
contributions to ejy, but because in some cases they can even allow some RH neutrino 
representations to have new decay modes that can change the way in which leptogenesis 
takes place. Once the conditions discussed above are satisfied not much more, from a 
general perspective, can be said and the way in which the B — L asymmetry is generated 
depends upon the particular flavor model. Hereafter are discussion will rely on a particular 
A4 flavor model realization. 

3 Setup 

We consider the non-supersymmetric version of a model inspired by the Altarelli-Feruglio 
model discussed in [TH] of which the type-I seesaw formulation has been analyzed in |2Uj . 
In the original model supersymmetry is introduced to induce the correct spontaneous 
breaking of the flavor symmetry. Here we assume that by adding additional discrete 
Abelian symmetries or ad hoc soft terms the scalar potential may be arranged in such a 
way that the desired breaking is realized. At this level the model presented is still a toy 
model, however our findings will hold even in its supersymmetric version. 

3.1 The A4 group 

Before entering into the details of the model, for completeness, we will briefly discuss the 
basic ingredients of the A4 discrete group in which the model presented here is based. A4 
is the group of even permutations of 4 objects. It has 4 irreducible representations: one 
triplet and three singlets 1, 1', 1". may be thought as generated by two elements S,T 
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satisfying 

S 2 = T 3 = {ST) 3 = 1 . (1) 

In what follows we will work in the At basis in which the triplet representation of T is 
diagonal, namely 

/ -1 2 2 \ / 1 \ 



2-12 
V 2 2 -1 J 



T 



u 
\0 cu 2 J 



(2) 



with cu 3 = 1. The multiplication rules in this basis are given by 

(afe)! = (aifti + a 2 b 3 + a 3 b 2 ) , 
(ab)v = (a 3 6 3 + aib 2 + a 2 h) , 
(a6)i« = (a 2 b 2 + ai6 3 + a 3 &i) , 

(ab) 3s = (2ai6i - a 2 ^3 - a 3 6 2 , 2a 3 6 3 - a 2 ^i - ai& 2 , 2a 2 6 2 - a 3 &i - ai& 3 ) , 
( a ^)3 a = (a 2 & 3 - a 3 6 2 , ai6 2 - a 2 &i, 03^1 - ai6 3 ) , 



(3) 



where a and 6 are triplets of A 4 , namely a ~ (ai, a 2 , a 3 ), 6 ~ (61, 6 2 , 63) and for the singlet 
representations the multiplication rules are trivial 



i' <g> 1" = 1 , v <g> 1 = i' , r (8i 1 = 1" . 



(4) 



3.2 The model 

In our model four RH neutrinos are added to the SM field content. Three of them, vt, 
form an A 4 triplet while the fourth, z/ 4 , is an A 4 singlet. The SM lepton doublets, l 1: l 2) / 3 , 
transform as an A A triplet. For simplicity we assign flavor charges using the Weyl spinor 
notation. The RH 4-dim Majorana fermion will be therefore defined as 



N = 




(5) 



It is important to notice that if I ~ (h,h,h) transforms as a triplet the requirement 
of recovering the correct kinetic term according to the group multiplication rules imposes 
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that ft transforms as a triplet but ordered as ft ~ 4, 4)- O n the contrary we order 
as (u\, ul, ul) thus v T ~ (z/i, z/ 3 , z/ 2 ). 

RH charged leptons transform as the 3 one dimensional representation of A4, namely 
1,1', 1". Two A 4 scalar triplets, <pT and <fis, our flavons, are added. Once the flavor 
symmetry is broken they will give rise to the correct mass matrices. 

3.2.1 Neutrino mass matrices 

Given the field content previously described the Lagrangian for the lepton sector — not 
including the kinetic terms — reads as 

-C = -t^{vtr v tr + VtrVtr)i + -i^Wr^r + Kr^ar)! 

+ A [^trVtr^s] 1 + A* [vtrVtr(P*s] 1 

+ t ["trM iKr + C v% [vtr4>s\ 1 

+ yi e a ^ TR lf) x H^ + yle^UTR^ 

+ yJ^wa^ + H-c. (6) 

We have assumed the presence of an Abelian with N > 2 that forbids <ps,T and 4>* ST 
to have the same couplings and we are assuming the flavons to be complex fields. Note 
that in this class of models this kind of Abelian symmetries are always present to prevent 
interferences between and (f>s- 

Note also that we can rephase the RH neutrinos to make M v t and M„ 4 real. After 
doing so we still have the freedom to rephase the lepton A4 triplet and the flavon A4 
singlet, absorbing in that way other two CP phases. These basis rotations allow us to 
choose A and y 3 to be real. 

In eq. ^ the last row describes charged leptons: i stays for 1,1', 1" and the [...], 
stays for the triplet contractions in the one dimensional representations. In the other 
rows [...]i stays for 2 or 3 triplets contracted in a singlet. Greek indices a and /3 label 
SU{2) L degrees of freedom, H the SM Higgs doublet, and as usual H = ia 2 H. We assume 
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that the additional symmetry forbids the coupling of 4>t and 05 to neutrinos and 
charged leptons respectively. The A4 basis chosen is useful because when <pT develops a 
vev according to (0^) ~ 1^(1,0, 0) the charged lepton Yukawa mass matrix is diagonal. 
On the other hand when <p s acquires the vev {(f> s ) ~ v s (l, 1, 1) the Dirac Yukawa matrix, 
Y U) and the RH neutrino mass matrix, Mn are diagonalized by the so-called TBM mixing 
matrix. Thus after electroweak symmetry breaking the light neutrino mass matrix, 



m,. 



-rrlF D ■ M N X -m D 



(7) 



is 



diagonalized by the TBM mixing matrix as well. Clearly m D = Y v vh, with (H) = vh- 
Both Mn and Y v get a contribution above (>) and below (<) the scale vf ~ vt ~ vs, 
so we may write 



M N = M N + M N , 



with 



M > 



N 



( A \ 
A 
A 



Y > _|_ y< 



(8) 



N 



\ D J 

( IB —B —B C \ 

-B 2B -B C 
-B -B 2B C 

c c c 



a o\ 
0a 
a 

y j 

( 2b -b 

-b 2b 

-b -b 

y c c 



V 



-b\ 

-b 



2b 

c 



(9) 



/ 



with A ~ M u t, B ~ Xvs, C ~ £i>s, -D ~ and a ~ yi, 6 ~ ^fs/A, c ~ y^vs/A for the 
RH neutrino and Yukawa Dirac mass matrix respectively. Without loss of generality v s 
may be taken real. 

Defining = fs/A, clearly in the limit — >■ the symmetry is restored and the light 
neutrinos are degenerate. Thus we expect in the majority of the cases a quasi-degenerate 
(QD) neutrino mass spectrum in which \yi\ 2 /M vT controls the absolute mass scale while 
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present Qvlfi bounds 




mi(eV) 



Figure 1: The predictions for \mpp\ as a function of the lightest neutrino mass. The 
natural spectrum predicted by the model is QD, with both normal ordering (NO) and inverse 
ordering (10) as indicated by the analytical approximations. Both QD-NO (QD-IO) and 
NH (IH) spectrum are indicated with green points (blue crosses). See the text for more 
details. 

1/2, A, es = vs/A parametrize the neutrino atmospheric mass splitting. We may find an 
approximate analytical solution for the spectrum expanding in es- We have 



v 

m ~ 1 2/i 



2 

2 V H 



~ 2 + 2 (£ * + £ *' <■■•><■■" o 2 > • («» 



(\yi 


2 i 2 -2\yx\y& cos <; 


6 W1 +y| cos20 w ) 




2/i I ( 


A-2 


2/2 


cos A0 12 ) 



with e x = (M v x - M u4 )/M uT , 4> yi 2 = Arg(j/i )2 ) and A0 i2 



'2/1 ry2 • 



Equation (10) holds only in the regime < ex- When ex < £5 the analytical expres- 
sions become more cumbersome because the QD scenario is broken and both, the normal 
hierarchical (NH) and inverse hierarchical (IH) neutrino mass spectra become possible. 
The spectrum predicted by our model is shown in figjljby means of mpp, the parameter rele- 
vant for neutrinoless double beta decay defined as mpp = [£/tbm diag(m 1; m 2 , m 3 ) U^ BM ]u- 
It has been obtained by numerically diagonalizing the neutrino mass matrix below the 
scale Vs, fixing the neutrino mixing angles according to the TBM scheme and requiring 
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the solar and atmospheric mass splittings to lie within their currents 3 a ranges [17J. The 
parameters entering in the mass matrix were varied according to: [10~ 3 , 1CT 1 ] for the di- 
mensionless parameters, es constrained to be below 10 _1 and RH neutrino masses in the 
range [10 7 , 10 14 ] GeV. For all the points A was taken to be larger than the heaviest RH 
neutrino mass. For completeness in the figure we have shown the future experimental 
bounds on \mpp\ and mj. 

For what concerns the flavon sector the Zn symmetries allows only the mass term 

M}(Msh ■ (11) 

Recalling now that (p* s ~ (4>*g 1 , <ft* S3 , 4>* S2 ) we have that the flavon mass matrix is diagonal 
and CP even and odd states are degenerate. 

3.2.2 Flavon Interaction Matrices 

Above Vp instead of .M^, Y"^ we have flavon-neutrino and flavon-Higgs- neutrino interac- 
tion matrices In,u- Starting from the interactions 

1 1 

+ ^ 2 e Q/3 [4^05]i^ + ^^e^[4 Q ^0s]i^, (12) 
we may write them using the 4-component spinors N and PlI = Il as 

+Z^j0sJ i i6^# i P i ^+Z^(^J^e a ^fP fl iV,#^, (13) 
where = Zjy ,Zq = Zq' and k labels the 0s, <f>* s flavons. Notice that in eq. (13) and 



for the rest of the paper we will indicate with Ni the four RH neutrinos of the model under 
study, while in sec. 2 Ni was referred to the group representations. 
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First of all we change basis going in the basis in which is diagonal. Thus we have 
M> = Ul-M>- U R = Diag(M„ T , M uT , M uT , M„ 4 ) = Diag (M N , 2 , 3) , M Na ) , 



U, 



R 



Y> 

1 D 



1 



JJ ] ■ Y> 



I V2 \ 
1 -i 
1 i 

o o o v 7 ^ y 

/v^ o\ 

1 1 

-i i 

\ / 



2V2 



(14) 



It proves useful to write the interaction matrices as 



z& a (& fc )« = (z&JS& 



and similarly for ^ L (<Ps k )ijj^'r> L ( < Ps k )ij- The X^, are 3x4 matrices, being 4 the total 
number of RH neutrinos. Equation ( 15 ) may appear as redundant since in our model 4>s 
couples always linearly. However the advantage of our notation is that it holds even when 
operators of dimension higher than 5 are included. 



3.3 CP asymmetries 

In the standard leptogenesis scenario the lightest RH neutrino CP asymmetry, €n, arises 
from the interference between the tree-level decay Feynman diagram and the one-loop ver- 
tex and wave function corrections [2TJ • Since N4 does not have renormalizable couplings to 
the lepton doublets such diagrams do not exist in the case under consideration, regardless 
of the RH neutrino mass spectrum. New contributions due to the presence of the flavons 
degrees of freedom exist and depend upon the RH neutrino spectrum: 
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(a) (b) 
Figure 2: Tree-level and one-loop correction diagrams accounting for e^. 

1. The > Mn 4 case: Ni has standard LiH tree- level decays and, given the inter- 
actions in the Lagrangian ([6]), the only possible correction to this process arise from 
the one-loop correction to the effective vertex N^iLiH, as shown in fig. [2] Thus, 
the CP asymmetry in this case is obtained from the interference between diagrams 
ga) andgb). 

2. The Mjv 4 > Mjy 4 case: Since N4 couples to N<pi at the renormalizable level a CP 
asymmetry for the two body decay process N± — » Nifa can be calculated from the 
interference between the corresponding tree-level diagram and the two-loop level 
diagram involving both effective couplings N^L^iH and NiL^iH (since iV 4 does 
not couples to lepton doublets at the renormalizable level the one-loop correction 
to the process N± — > Ni<f>i does not exist). There is another option involving N4 
three-body decays induced by the effective coupling N^LifaH. In this scenario a 
one-loop correction to the effective process does not exist either and the calculation 
of the CP asymmetry relies again on the two-loop level correction of the previous 
case. 

In case[T]the CP asymmetry arises in a different way compared to the standard case but 
in what regards the generation of the B — L asymmetry there is no difference. In contrast, 
the cases in [2] are quite different: for the three-body decay scenario the differences are 
obvious, for the other scenario leptogenesis will take place in two stages, a first stage in 
which an asymmetry in iVj is generated via the decays N 4 — > and a second stage in 
which the asymmetry in Ni is partially transfered to the lepton doublets via iVj decays 
and scatterings (a scenario of this kind has been discussed in P, [22]). Note that in this 
case the CP asymmetry, being a two- loop order effect, would most likely yield a very tiny 
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B — L asymmetry. All these scenarios however exhibit a common feature, the generation 
of a B — L asymmetry takes place in the flavor symmetric phase. So from now on we will 
focus on case [TJ that as was already pointed out resembles standard leptogenesis. 
The CP asymmetry in the decay of Aj is defined according to 

k=e,/i,T k=e,fj,,r 

where r^- fe (T^ h ) denotes the A; partial decay width for final states of flavor k and carrying 
+1 (—1) unit of lepton number, and e^ h are the flavored CP asymmetries. However, here 
we are working in the context of an exact flavor symmetry. Since flavor is unbroken only 
flavor conserving processes may happen, that in our framework means k = i. Moreover 
flavor invariance and the representation used for our right and left-handed neutrinos imply 
that the three RH neutrinos produce the same amount of CP asymmetry and have exactly 
the same dynamics. Due to the complex nature of the scalar field components fa running 
in the loop |2|b) there is only one possible one-loop diagram of that type (contrary to the 
standard leptogenesis case for the wave-function correction), so the interference between 
the tree and one- loop level amplitudes (A^o and A^i) involves only one term. For two- 
body decays this interference is phase-space independent and consequently the calculation 
of ej^ can be simply done in terms of the products of Aio and Ai± and approximating the 



denominator of (16) with |7W 1 2 In the limit M N . » M Ni ,M^ the CP asymmetry 
can be written as 

<_>/< | -i ± 

In the general case, without assuming any large hierarchy among the heaviest RH neutrino 
representation, the lightest one and the flavons, the expression for is far more involved. 
Fig. [3j obtained from the exact expression, shows the possible values of the CP asymmetry 
as a function of the effective cut-off scale A. 



3.4 Generation of the B — L asymmetry 

The generation of the B — L asymmetry is entirely determined by A, dynamics but its final 
value depends on the washout induced by the At flavor singlet, A 4 . Thus, leptogenesis in 
the case we are interested in is a two-step process: generation of the B — L asymmetry 
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10 10 10 10 11 10 12 10 13 10 14 10 15 
A (GeV) 

Figure 3: as a function of the cutoff scale A for the spectra shown in fig. [I| Green points 
(blue crosses) corresponds to QD-NO and NH (QD-IO and IH) neutrino mass spectra. 
The input parameters are generated as in fig^ is computed according to eq. (16) and 
eq. (F 



and its subsequent washout via A/4 interactions (such scenario has been analysed in the 
context of type-Ill seesaw in [23] ). In what follows we will analyze both stages in the 
unflavored regimen. 



3.4.1 Ni dynamics 

The determination of the B — L asymmetry relies on the solution of the kinetic equations 
for the Ni abundance and the B — L asymmetry itself. At leading order in the coupling y 1; 
that is to say including only Ni — > LH decays and inverse decays and AL = 2 scatterings 
(LH^ -H- LW and LH -H- LH) PI they can be written according to 



dY, 



N, 



1 



dzi 



dYl 



dzi 



s(zi)H(zi)Zi 
1 

s(zi)H(zi)zi 



Yn, 



1 7a(*)» 



1 \e Ni + 



Y A 



B-L 



Eq 



(19) 



Y^izi) ) - ' 2Y 

4 The inclusion of these processes is mandatory to obtain kinetic equations with the correct thermody- 
namical behavior [5S1 [2H ■ 
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where Zi = M N ./T, s is the entropy density, Y x = n x /s (with n x the X number density), 
Yl = 2Ye + Y e (the lepton asymmetry distributed in left-handed and RH degrees of free- 
dom), H(z) is the expansion rate of the Universe and the reaction density 7Dj(^) is given 
by 

. 1 M% Ki(zj) _ 

7a (*) = ^3 ^ > ( 2 °) 



with v ~ 174 GeV, i^i(^) the modified Bessel function of first-type and the parameter 



ttlt = v 2 \yi\ 2 1 'M^. An exact solution of the kinetic equations in (18) and (19) can only 
be done numerically, however a reliable approximate solution can be found (26] , which we 



now discuss in turn. Equations (18) and (19) can be recasted according to 

dY Ni 



dzi 



art 



dzi 



-D T { Zi 
-e Nz D T (z 



Y^, [Zi 



YnM) 



W? D ( Zi )Y Al 



where the new decay and inverse-decay functions read 
D T (zi) 



K -r^irr\ and W f D (Zi) = - A K T zlK 1 (z i ) 

J<2{Zi) 4 



(21) 



(22) 



with K T = m T /m* (m* = 8irv 2 H(zi = l)/M%. = 1.08 x 10~ 3 eV) and K 2 (Zi) is the 
modified Bessel function of the second-type. In terms of Kt the strong (weak) washout 
regimen is defined as Kt ^> 1 (Kt < 1). 



The B — L asymmetry is obtained from the formal integration of eqs. (21) by means 
of the integrating factor technique: 



Yl B _ L (^) = -3 x e Ni Y^( Zl -> 0) fa) 



Here r)(zi) is the efficiency function defined as 



1 



-> 0) 



' d Z ' dY NX z ') r - Q dz"Wf D {z") 



2() 



dz' 



(23) 



(24) 



Note that we have included a factor of 3 in (23) to account for the N{ flavor degrees of 



freedom. The final B — L asymmetry is therefore obtained for z\ — > oo once the parameters 
K T and e Xi are specified. 

The problem of determining the final Y£ analytically is thus reduced to find an 
approximate expression for the efficiency function at Zi — > oo (efficiency factor). Such 
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Figure 4: B — L asymmetry produced by the flavor triplet (N{) dynamics. The values 
€t = 10~ 6 — 10 _1 correspond to extreme cases. 

an expression can be derived in the strong washout regimen by: (i) noting that at 
low temperatures Y N .(zi) follows closely the equilibrium distribution, so the replacement 



dY N .(zi) / dzi — > dY^(zi) / dzi in eq. (24) can be done; (ii) replacing the washout function 
^w{ z i) by Win(zi) = zb Wf D {zi)J z, where zb is the minimum of the function 



In 



W 



ID\ 



+ 



dz"Wj D {z"). 



(25) 



Following this procedure the efficiency factor can be derived 

2 



V 



K T z B (K T ] 



(i 



-K T z B (K T )l2\ 



with 



zb(K : 



r 



bn ^ 
2 I 1024 



In 



3125 ttK 



T 



1024 



(26) 



(27) 



With eqs. (26) and (27) at hand we can determine the maximum and minimum (still 



consistent with the measured baryon asymmetry) B — L asymmetry one can get through 
Ni dynamics. The results are displayed in figure |4} Particularly relevant is the maximum 
value Y£ ~ 10~ 5 as it allows to derive an upper bound on the washout induced by 
dynamics. 
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L 



~- H 



Figure 5: Relevant 1 f> 3 and 2 f> 2 s, t and u scattering processes accounting for the A 4 
singlet washouts. 

3.4.2 N 4 washout 

The B — L asymmetry produced at Zi ~ 1 remains frozen up to the temperature at which 
iV 4 washouts become effective, z 4 = M Ni /T ~ 1. Since iV 4 couples to lepton doublets 
via an effective five-dimensional operator the dynamics of N± washouts, at leading order 
in the couplings, involves not only the processes iV 4 -h- L<pH but the 2 -h- 2 s, t and u 
channel scatterings (see figure [5]), in contrast to the standard leptogenesis scenario. The 
derivation of the corresponding kinetic equations in this case is tricky and requires -even 
at leading order in the couplings- the inclusion of 3 -H- 3 and 2 -H- 4 scattering processes 
(see ref. [10J for more details). Since ejv 4 = the kinetic equations accounting for the iV 4 
washouts can be written according to 



where tI^C 2 ^) ^ s the reaction density for the 2 -H- 2 s-channel scattering process and 
TtatO&i) involves the reaction densities for the full set of processes shown in figure [5l namely 




(28) 



(29) 




(30) 



C=s,t,u 



As explained in appendix [X] all the reaction densities can be written in terms of the total 
decay width r(A^ 4 — > LHtfi), which we have calculated to be 



r(iV 4 -> LH(j>) 





(31) 



192tt 3 A 2 
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10" 



10" 



10 



10" 



10 



Figure 6: Reaction densities as a function of Zi for the different processes present in 
washout (left panel) and total reaction densities for different values of the decay parameter 
Ks (right panel). As in the standard case the strong washout (weak washout) regimen is 
defined according to K$ ^> 1 (Ks <lj. 



In terms of the reaction densities given in (43 ) the kinetic equations in (28 ) can be rewritten 



in such a way they resemble eqs. (21 ): 

dY Nt 



dzA 



-D s (z 4 ) 



Y Ni (zi) 



dz^ 



-W? D {z A )Yl B _ L , 



(32) 



where now the functions D s and Wf D are given by 



K s z\ 



Wf D {z i ) = -K s zl 



K^ + ^S^ + St^)) 



K 1 {z 4 ) + 



3 / Y A 



■ N 4 



(33) 



Some words are in order regarding these equations. The relativistic degrees of freedom 
are = 118, as in our calculations we use Maxwell-Boltzmann distributions, and the 



functions S s t are given in eqs. (44) in the appendix. The decay parameter Ks is defined 



in the same way it is defined in the case of Ni dynamics, Ks = ms/m+ but with 



m s 



Q7TV 



r(iV 4 -> LH(j)) . 



(34) 




10" 

M Ni [GeV] 



Figure 7: Washout induced by N4 on the maximum Y£ (see fig. 4) generated in Nt 
dynamics (left panel). The largest allowed K$ for which the resulting Y& B still fits the 
measured value is K^ iax = 4.6 x 1CT 2 , any value for which K$ > K™ ax is excluded. On the 
right panel allowed regions of A — M^ s as required by the condition K$ < K™ ax . 



The presence of the 2 — > 2 scattering processes may drive the system to the strong washout 
regimen even when the 1 — > 3 process is slow. Thus, the appropiate definition of the strong 
(weak) washout regimen in this case reads: 

7tot0i) 



s{Zi)H(zi)Zi 



> 1 (< 1 for weak washout) . 



(35) 



Figure [6] (left panel) shows an example in which though "fx-^^Zi) J s(zi)H(zi)Zi\ Zi ~\ <^ 1 
the system is driven to the strong washout regimen by scattering processes. Note however 
that the condition K$ 3> 1 (K$ <C 1) still determines the regimen in which the washout 
dynamics of takes place, as can be seen in figure [6] (right panel). 

can be 
The 



From the integration of eqs. (32) an upper bound on K$ for a given Y£ 



B-L 

-10 



determined by the condition of not erasing this asymmetry below ~ 2.6 x 10" 
maximum value i^™ ax is found for the largest possible B — L asymmetry generated in TV, 
dynamics, that as has been argued in sec. 3.4.1 we have found to be ~ 10~ 5 . Figure [7] (left 
panel) shows the final Y£ b _ l matches the required value ~ 2.6 x 1CT 10 (for Y£ = 10~ 5 ) 
when i^™ ax ~ 4.6 x 10 -2 , any value K$ > K™ ax will induce a washout that will damp the 
B — L asymmetry below the allowed value. 
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Taking \y 3 \ = 10 2 the decay parameter K s becomes 

— 9 (^)(T) 2 < 

with the purpose of placing the more stringent bounds on the A — Mn a plane we fix 
10 2 < A/Mjv 4 < 10 4 and take into account the restriction Kg < -K]? ax . The results are 
displayed in figure [7] (right panel) where the allowed A — M Ni region can be seen. Any 
discussion of leptogenesis in the scenario we have considered here should be done at least 
within that region. 



4 Conclusions 

In this paper we have study the necessary conditions that have to be satisfied whenever the 
generation of the cosmic baryon asymmetry of the Universe via leptogenesis takes place in 
the presence of a lepton flavor symmetry accounting for lepton mixing. In the scenario we 
have discussed, leptogenesis occurs in the flavor symmetric regime thus before the flavons 
(that trigger the breaking of the flavor symmetry) acquire vevs, accordingly the decays 
responsible for generating a net B — L asymmetry are liable of selection rules dictated by 
the flavor symmetry. 

In the core of the paper we exemplify how the general conditions for the generation of 
the baryon asymmetry, in the flavor symmetric phase, work by analysing a specific model 
based on the flavor symmetry A±. We briefly discussed the low energy phenomenology 
of the model and studied in detail, by using the corresponding kinetics equations, the 
generation of the baryon asymmetry. In the model considered, due to the constraints 
imposed by A4, the asymmetry proceeds through the CP violating and out-of-equilibrium 
decays of the heaviest RH neutrino A 4 representation. Subsequent washouts induced 
by the lightest A4 representation, being potentially dangerous, were properly taken into 
account. Our onset shows these washouts can always be circumvented and the correct 
amount of baryon asymmetry can be produced. 

In conclusion we have shown that under certain conditions, in models containing fla- 
vor symmetries in the lepton sector, leptogenesis can occur even in the flavor symmetric 
phase. The conditions we have enumerated can be regarded as a general recipe for con- 
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structing lepton flavor models in which the lepton number violating scale is above the 
flavor symmetry breaking scale and the generation of the baryon asymmetry proceeds via 
leptogenesis. 
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A Appendix: Reaction densities for 1 — >• 3 and 2 — > 2 
processes 

In this appendix we present the relevant equations used in the calculations discussed in 



section 3.4.2 The thermally averaged reaction densities for 1 — > 3 and 2—7-2 processes 



are given by [25] 

Eq Ki(z4) 

7i^3 = n N H ri_> 3 , (37) 

t£> 2 = g , q 4 / dxy/xK 1 (z i y/x)a c (x) } (38) 

where x = s/M^ 4 (with s the center of mass energy) C = s,t,u, r\^ 3 = T(A^ 4 — > LH<p) 
and a(x), the reduced cross section, defined as 

d{x) = 2M^ 4 xA(l,x" 1 ,0)a(i) with A(a, b, c) = (a - b - cf - 46c . (39) 

Neglecting the lepton doublets, Higgs and flavones masses we have found for the differential 
cross sections the following results: 

da s _ 1 |y 3 | 2 1 1 
~dT ~ 16tt "A^" M* 4 1-x ' 
da 1 1 \y,\ 2 lift 



dt 16vr A2 M^(l-x) 2 V M 2 Na 

d^ = ^\y£^ l —(x + ^-\ (10) 
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Integrating over t in the range t_ = (1 — x) and t + = and using the definition for 



the reduced cross section, eq. (39), we get 



irf (m n ,v [x-if 



a J * • (41) 
s ..« w= M(^y(^. ,42) 

v ; 16tt V A J x K J 

With these results at hand and taking into account the expression for Ti^ 3 given in eq. 



(31) the different reaction densities become 



i V4 




r 


7T 2 ' 


z 




S s (z) 




2tt 2 


z 




3M^ 4 







4tt 2 

where the functions S s ,t,u{z) are given by 



71^3 

7^2 - —^T — 1 L-:{ 

7 2 U 2 = -7-^ — — r i->3 , (43) 



= / dx K - — —>-K x {zyfx) and = / dx—^ K x {z>/x) . (44) 
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